Magnetic Correlations at Graphene Edges 
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Magnetic zigzag edges of graphene are considered as a basis for novel spintronics devices despite 
the fact that no true long-range magnetic order is possible in one dimension. We study the transverse 
and longitudinal fluctuations of magnetic moments at zigzag edges of graphene from first principles. 
We find a high value for the spin wave stiffness D — 2100 meV A 2 and a spin-collinear domain wall 
creation energy iJdw = 114 meV accompanied by low magnetic anisotropy. Above the crossover 
temperature T x ~10 K the spin correlation length £ cx T -1 limits the long-range magnetic order to 
~1 nm at 300 K while below T x it grows exponentially with decreasing temperature. We discuss 
possible ways of increasing the range of magnetic order and effects of edge roughness on it. 

PACS numbers: 75.75.+a, 75.40.Cx, 81.05.Uw, 85.75.-d 



Graphene, a two-dimensional form of carbon, has at- 
tracted considerable attention due to its unique phys- 
ical properties and potential technological applications 
0,0. The possibility of desi gning graphene-based mag- 
netic nanostructures is particularly intriguing and has 
been fuelled by the recent experimental observations of 
magnetism in graphitic materials 0, 0). A number of 
exceptional nanoscale spintronics devices built around 
the phenomenon of spin polarization localized at one- 
dimensional (ID) zigzag edges of graphene have been 
proposed d, 0, U However, feasibility of such de- 
vices is questioned by the fact that no true long-range 
magnetic ordering in ID systems is possible at finite tem- 
peratures 0. Nevertheless, nanometer range spin corre- 
lation lengths in certain ID systems have been achieved 
in practice [loj | . Establishing the range of magnetic or- 
der at graphene edges as well as the underlying physical 
mechanisms is thus crucial for practical realization of the 
proposed spintronics devices. 

In this Letter we study the magnetic correlations at 
zigzag edges of graphene by investigating the transverse 
and longitudinal fluctuations of magnetic moments from 
first principles. While the transverse excitations (spin 
waves) are characterized by the continuous rotation of the 
electron spin moments along the edge (Fig. QJ,), the lon- 
gitudinal fluctuations affect the spin correlation length 
only if an inversion of magnetic moments resulting in ap- 
pearance of a spin-collinear domain wall [H[ takes place 
(Fig. |TJd) . The evaluated energies of these low-energy 
excitations mapped onto the classical Heisenberg/Ising 
models allow us to estimate the spin correlation lengths 
at different temperatures. Finally, possible ways of in- 
creasing the spin correlation length and the effects of edge 
roughness are discussed. 

The first-principles calculations of the magnetic ex- 




FIG. 1: (Color online) Schematic representation of the trans- 
verse (a) and longitudinal (b) low-energy spin excitation at 
graphene zigzag edges. The magnetic moments of the out- 
ermost edge atoms are shown by arrows. The direction of 
magnetic moments is represented by direction and color of 
the arrows while the magnitude is illustrated through the ar- 
row lengths and color intensities. 



citations are performed on the density functional the- 
ory (DFT) level using the Perdew-Burke-Ernzcrhof 
exchange-correlation functional [l2T ] . A non-collinear spin 
DFT formalism [H, [T3] implemented in the PWSCF plane 
wave pseudopotential code 11511 in combination with the 
ultrasoft pseudopotentials [1 61 ] and a plane wave kinetic 
energy cutoff of 25 Ry is used to study spin wave modes. 
Much larger supercells are required to obtain converged 
results for the spin-collinear domain walls. These calcu- 
lations are performed using the standard spin-polarized 
DFT scheme implemented in the SIESTA code [ijj to- 
gether with a double-^ plus polarization basis set, an en- 
ergy cutoff of 200 Ry and normconserving pseudopoten- 
tials [HI . Test calculations performed on limited size sys- 
tems verify that both codes provide results in close agree- 
ment. The model systems considered are the hydrogen- 
terminated periodic one-dimensional graphene nanorib- 
bons of different widths and supercell lengths relaxed in 
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FIG. 2: (Color online) (a): Spin density isosurface plot for the 
collinear domain wall excitation at a zigzag edge of graphene. 
Spin populations m (b) and spin-resolved projected density 
of states (c) for the outermost edge atoms. The projected 
density of states values for spin-up and spin-down electrons 
are indicated by the intensities of red and blue colors, respec- 
tively. The edge atoms are numbered with n. 



their ground state configurations. 

The ground state electronic configurations of zigzag 
graphene nanoribbons is characterized by the ferromag- 
netic arrangement of spins along the edges and antifer- 
romagnetic coupling of the spins at the opposite edges 
[l9l | . To obtain a spin- wave-excited state we perform 
constrained self-consistent calculations with a penalty 
functional term [2(| added to the total energy expression 
in order to induce small non-collinear deviations of the 
magnetization directions from the spin-collinear ground- 
state configuration. The total energy energy difference 
is mapped onto the quadratic spin-wave dispersion re- 
lation, E(q) = nq 2 , with k = 320 meV A 2 . At a zigzag 
edge of graphene the magnetic moments of the outermost 
edge atoms m c d gc = 0.28 /ib while the magnetic moments 
localized on the atoms belonging to the A and B sublat- 
tices within a single edge unit cell are to a = 0.43 and 
tob = —0.13 fiB, respectively. This yields a total mag- 
netic moment of to = toa + tob = 0.30 ^b per unit cell of 
zigzag edge. The obtained value of to agrees with the fact 
that in zigzag graphene nanoribbons a flat band develops 
in one-third of the ID Brillouin zone (27r/3 < \ka z \ < tt; 
a z = 2.46 A is the unit cell length) when electron-electron 
interactions are not taken into account [2l|. The spin- 
wave stiffness constant D — 2n/m turns out to be 2100 
meV A 2 . Actually, this is a very high value which is 
about one order of magnitude higher than the stiffness 



Thus, our results confirm the expectation of higher spin 
stiffness values in magnetic materials based on sp ele- 
ments compared to d element materials 24] . 

In sp-electron itinerant-electron magnets, Stoner-type 
longitudinal spin fluctuations may be essential 24j. To 
estimate their characteristic energy we study collinear 
domain walls at the graphene zigzag edge. We have per- 
formed the calculations on a large graphene nanoribbon 
supercells (up to ml. 8 nm wide and 6 nm long). In or- 
der to converge the self-consistent calculations to the do- 
main wall solution we provide an appropriate initial mag- 
netizations of edge atoms with two equidistant domain 
walls per unit cell for maintaining periodicity along the 
nanoribbon direction. Figure [5^ illustrates the distribu- 
tion of the spin density at such a domain wall located in 
the center of the edge fragment shown. The spin popu- 
lations of the outermost edge atoms (Fig. [2Jd) show that 
the domain wall is practically localized within two unit 
cells (0.5 nm) and the magnetization exhibits weak os- 
cillations close to the kink. The spin-resolved projected 
density of states for the outermost edge atoms (Fig. (2J:) 
shows an avoided crossing pattern with band gap dimin- 
ishing (but not closing) at the domain wall. From the 
total energy difference we find a collinear domain wall 
creation energy Edw = 114 meV per edge. 

In order to determine the magnetic correlation param- 
eters in the presence of spin wave fluctuations we recall 
the nearest-neighbor ID classical Heisenberg model 
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where is the magnetic moment unit vector at site i and 
H is the external magnetic field vector. The Heisenberg 
coupling a — 2«/a 2 = 105 meV corresponds to the value 
of k calculated above from first principles. The axial 
anisotropy parameter d is expected to be small due to in- 
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trinsically weak spin-orbit coupling in graphene 
We obtain an-order-of-magnitude estimate for the mag- 
netic anisotropy d/a — 10~ 4 using the spin-orbit coupling 
strength of ~0.01 meV 25J predicted for graphene with 

The 



weak corrugations observed experimentally [271 . 
estimated d/a agrees with the recent measurements of 
2D magnetic correlations in irradiated graphite and 
with the electron spin resonance g-tensor anisotropies in 
molecular graphitic radicals 

Sim. 

The spin correlation length (a = x, y, z) de- 
fines the decay law of the spin correlation function 
(sfsf +l ) = (sf sf)exp(—l I £ Q ) , i.e. the range of magnetic 



constant of bec iron 22, k3J|, a three-dimensional solid 



with much larger magnetic moment of 2.2 /zb per atom. 



order. First, we evaluate the zero-field spin correlation 
length due to the transverse spin fluctuations as a func- 
tion of temperature (see Fig. [3]) 1321. Above the crossover 
temperature T x = \fad ~ 10 K [33| the small anisotropy 
term of the model Hamiltonian has practically no in- 
fluence and the system exhibits behavior typical for an 
isotropic Heisengerg model [36[ with ^ w w 300/T [nm] 
and (sfsf) = 1/3. Below T x the anisotropy term starts 
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FIG. 3: Correlation lengths of magnetization vector compo- 
nents orthogonal (£ 2 ) and parallel (£ x , £ y ) to the graphene 
plane as a function of temperature T for weakly anisotropic 
(d/a = 10 -4 ) and isotropic (d/a = 0) Heisenberg models. 



playing an important role and the solution exhibits a 
characteristic for ID Ising model exponential divergence 
of 4 2 w oc exp(V8ad/fcr) and (sfsf) = 1 for T -> K. 
The spin correlation length at zero field in the pres- 
ence of spin-collinear domain walls is the one for ID 
Ising model, £j* w w exp(E dw /kT). Since E dw > V8ad 
[3^ | the overall spin correlation length £ in the pres- 
ence of both transverse and longitudinal fluctuations, 
= CJ + Cdw ~ CJ is defined predominantly by the 
spin wave disorder. 

At room temperature (~300 K) the spin correlation 
length £ = 3.7 unit cells (~1 nm). This result implies 
that a spintronics device based on magnetic graphene 
edges can be operated at room temperature only if its di- 
mensions do not exceed several spin correlation lengths, 
i.e. several nanometers. The device dimensions can be 
scaled linearly by lowering the operation temperature 
and below T x this size could be extended beyond the mi- 
crometer scale. These estimations may first look rather 
disappointing, but nevertheless they are comparable to 
one of the most appealing example of ID magnetism: 
monoatomic Co chains on Pt substrate characterized by 
a ferromagnetic order range of «4 nm at 45 K [loj ]. In 
this d-element system ferromagnetic order stems mainly 
from the anomalously high magnetic anisotropy which is 
absent in graphene nanostructures. However, the lack of 
anisotropy is partially compensated by the high spin stiff- 
ness which results in considerable spin correlation lengths 
even in the isotropic regime above T x . While the spin 
stiffness constant can hardly be increased we suggest sev- 
eral ways of increasing the magnetic anisotropy (and thus 
T x ) by strengthening the spin-orbit coupling by increas- 
ing curvature, applying external electric field or coupling 
graphene to a substrate (H)]. Alternatively, the mag- 
netic anisotropics can be increased by chemical function- 
alization of graphene edges with heavy element functional 
groups (e.g. iodine) coupled to the spin-polarized edges 
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FIG. 4: Ideal zigzag edge of graphene (a) and various types of 
edge defects: missing or rehybridized edge atom (b), Stone- 
Wales defect (c), edge step (d), and 120° edge turn. The 
domain wall creation energies at these structures are shown. 



states via the exchange polarization [3J, [35j . Augment- 
ing the crossover temperature above 300 K would result 
in a significant increase of £ 2 to the length scales of the 
present-day semiconductor technology. 

Thus, the graphene edges at finite temperatures are 
not actually ferromagnetic but superparamagnetic ones. 
For the isotropic Heisenberg model the enhancement fac- 
tor for the susceptibility in comparison with one of non- 
interacting spins reads [36| 
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where u = coth(a/T) — T/a and the approximation being 
valid at a > T. At room temperature the susceptibility 
enhancement factor x/xo ~ 8. 

Although we found a relatively high value of E dv/ , the 
localized domain walls may become energetically more 
favorable at edge defects, and therefore we discuss cre- 
ation of localized domain walls at different types of topo- 
logical imperfections at zigzag graphene edge classified 
as shown in Fig. [4] The simplest case of edge roughness 
is a boundary atom missing from the 7r-conjugation net- 
work (Fig. [4j>) . Such sp 2 -vacancy formation may result 
from the rehybridization of an outermost atom into the 
sp 3 state due to chemical modification or because of the 
creation of a true vacancy. The domain wall creation en- 
ergy at an sp 3 -hybridized atom is found to be 24 meV, 
i.e. factor of 5 smaller than E dw = 114 meV for the ideal 
zigzag edge. Such decrease will have a dramatic effect 
on the long-range magnetic order at room temperature 
since -Edw is lowered to kT («25 meV at 300 K). An 
even more dramatic decrease to 4 meV is observed at the 
Stone- Wales defect (Fig. 2b) , a topological structure ob- 
tained by the 90°-rotation of a single C— C bond which 
locally breaks the bipartite lattice symmetry. The pres- 
ence of an edge step (Fig. [4b!) has a less severe effect and 
reduces E dw to 62 meV. A completely different situation 
is observed for a 120°-turn of the zigzag edge (Fig. 2b)- 
The antiferromagnetic arrangement of spins at the edge 
segments separated by the 120°-turn is by 22 meV more 
stable than the ferromagnetic arrangement. This is due 
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to the change of bipartite sublattice to which belong 
the outermost edge atoms and due to the antiferromag- 
netic coupling between the magnetic moments in different 
sublattices 37, HI]. Similar behavior has recently been 
pointed out for the edges of hexagonal graphene nanois- 
lands [39]. Domain walls are thus naturally pinned to 
such turns, although the energy difference is close to kT 
at room temperature. A "spin-inverter" device design 
based on such a 120°-turn topology can be anticipated. 
Simple chemical modifications which do not perturb the 
7r conjugation network at graphene edges show almost no 
effect on i?dw For an ideal zigzag edge terminated with 
electronegative fluorine atoms we find E^w = 117 meV 
very close to the value for the hydrogen-terminated edge 
(114 meV). 

To conclude, we have studied from first principles the 
energetics of transverse and longitudinal spin fluctuations 
at the one-dimensional magnetic zigzag edge of graphene. 
The transverse fluctuations characterized by the high 
spin stiffness constant are the main limiting factor of the 
spin correlation length which is found to be ~1 nm at 
room temperature. For the temperatures above ~10 K 
the spin correlation length is inversely proportional to the 
temperature due to the low magnetic anisotropy of the 
system. Below the crossover temperature the spin cor- 
relation length grows exponentially with decreasing tem- 
perature. We propose several approaches for extending 
the range of magnetic order by increasing the magnetic 
anisotropy in this carbon-based system and discuss the 
effect of edge roughness on the spin correlation length. 
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